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Abstract 


This  report  (Preface  PI,  P2;  Parts  A-E)  is  essentially  a 
Summary  of  the  five  M.S.  Theses  and  Ph.D.  Dissertations  completed 
during  the  contract  year  for  contract  No.  DA-49-146-XZ-494.  The 
titles  of  these  theses  and  dissertations  are  the  heading  titles 
for  each  of  the  five  Parts  A-E. 

This  Final  Report  is  more  than  a  mere  Summary  since  all  results 
have  been  condensed,  compared,  corrected  and  interpreted  as  a  group. 
Figures  A8,  C14  and  D4  exhibit  count-rates  and/or  particle  flux 
densities,  both  for  the  Natural  and  the  ARGUS  Events  1  and  2  radia¬ 
tions,  in  3-dimensional  B-L  coordinates.  These  3-dimensional  plots 
exhibit  most  clearly  the  coordination  of  all  the  data  and  are  pre¬ 
sented  for  the  first  time  in  this  Summary  Final  Report. 


PI.  Preface  -  Glossary 


SYMBOL 


DEFINITION 


A  Nucleon  number  (Ref.  PI,  Sec.  1) 

A 

A  Vector,  along  axis  of  minimum  moment  of  inertia  Ex¬ 

plorer  4.  Also  the  direction  of  the  satellite  symmetry 
axis.  (Ref.  Pl}  Fig.  2.2). 

Ae^£  Effective  area  (product  of  counter  efficiency  and  geg- 

4i  metric  area  A)  of  directional  counter  Explorer  4  (cm  ) 


a. 

l 


Scale  height  for  atmospheric  model  defined  by 

1  Ni 

In  r; -  ,  where 


a.  = 

l 


h.  h. 

l+l  l 


N 


i+1 


a 


A 

P 


A 

B 


A 


B 

o 


B 


m 


Bj. 


N.  *  number/unit  volume  at  altitude  h.  where  h.  is  the 
altitude  in  kilometers  (Ref.  P2)  1  1 

Equatorial  radius  of  international  ellipsoid  (km) 

Daily  equivalent  planetary  amplitude  (Ref.  P3,  Sec.  1.5) 

Three-hour  equivalent  planetary  amplitude  (Ref.  P3, 

Sec.  1.5) 

Geomagnetic  flux  density  vector  (gauss) 

Geomagnetic  flux  density  unit  vector,  B/B 

Magnetic  flux  density  unit  vector  at  the  time  the  satellite 
inertial  coordinate  system  is  established  (time  of  any  rf- 
null  prior  to  data  to  be  analyzed) 

Mirror  point  value  of  B  (gauss) 

A  vector  in  the  plane  perpendicular  to  the  magnetic  flux 
density  vector  B  (gauss) 


Bv,Bw,B^  Components  of  a  unit  vector  B  in  Vernal  Equinox  coordi- 

1  nates 

Be  Equatorial  value  of  B 

b  Polar  radius  cf  international  ellipsoid 

A 

b.,b0,b  Components  of  B  on  the  Xs,Ys,Zs  axes  0 f  satellite 

°  inertial  coordinate  systlm  (Ref.  PI,  Eq.  (3.14)) 


Pl-1 


c 

c 

c 

we> 

^obs 

C(6) 

c(t) 

E 

a 

E 

e 

f(h') 

F(oL,<f))J(l 

(f) 

H 

h>hr 

h 

I 

1(E) 


i(H) 

I  M 
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Normalized  (divided  by  B)  components  in  topcentric  co¬ 
ordinate  system  of  ephemcris  (north,  east,  down) 

Omnidirectional  count-rate  corrected  for  deadtime  (sec  j 

Speed  of  light  in  free  space 

Observed  count-rate  corrected  for  counter  dead-time  (sec"*) 

Observed  directional  count-rate  for  detector  A  (sec-*) 

Omnidirectional  observed  count-rate  (sec-1) 

True  directional  count-rate  (0^(0)  c^rregted  for  deadtime, 
efficiency  and  effective  area)  (see-  cm  ) 

Auto-covariance  function  (Ref.  PI,  P4,  Eq.  (A1.2)) 

Kinetic  energy  of  a  particle 

Electric  field  vector 

Electronic  charge 

Mirror-point  density  (Ref.  P2,  Eq.  (2.24)) 

Fraction  of  particles  with  directions  within  dCL  at 
(Ref.  P3,  Sec.  4.S) 

Fraction  of  fission  electrons  that  penetrate  a  shielding 
(Ref.  P3,  Sec.  2.4) 

Receiver  antenna  height  above  the  international  ellipsoid 
(mean  sea  level)  (earth  radii) 

Relativistic  Hamiltonian  (Ref.  P3,  Sec.  3.2) 

Distance  measured  from  the  center  of  the  earth  to  a 
field  line  (km) 

Parameter  related  to  the  longitudinal  invariant  (Ref.  P3, 

Eq.  (3.22)) 

Omnidirectional  intensity  of  particles  with  energy  greater 
than  E  (Ref.  P3,  Eq.  (4.38)) 

Omnidirectional  spectral  intensity  (Ref.  P3,  Eq.  (4.32)) 
Modified  Bessel  function  of  order  p 
Moment  of  inertia  about  the  X'  axis 
Moment  of  inertia  about  the  Y'  axis 


Pl-2 


Jj(x) 

j(c6,^>;0) 

J(O') 

J 

o 

i 
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m 

m 

o 

N 

N(E) 


n.  ,-n. 
l+l  1 


n 

j*. 

P 


/s 

P 


P(f) 


r 

R 


Moment  of  inertia  about  the  Z'  axis 

First  order  Bessel  function  regular  at  the  origins 

Unidirectional  spectral  intensity  (Ref.  P3,  Eq.  (4.32)) 

Directional  flux  density  of  particles  as  a  function  of 
angle  between  J  and  a  plane  perpendicular  to  B 

Omnidirectional  flux  density  (particles  cm~  sec"  ) 

Length  of  nose  cone  and  instrument  section  Explorer  4 

Angular  momentum  unit  vector 

Mcllwain's  geomagnetic  shell  parameter  (earth  radii) 

Earth  dipole  moment  (Ref.  P3,  Eq.  (1.9)) 

Relativistic  mass  of  electron 
Rest  mass  of  electron 

The  number  of  digitizer  counts  per  digitizing  interval  T 

Absolute  differential  spectral  intensity,  electrons/fission/MeV 
(Ref.  P3,  Eq.  (1.5)) 

The  number  of  digitizer  counts  measured  to  correspond  to 
one  frequency  change  (flip)  on  data  record 

Index  of  refraction  (Ref.  P3,  Eq.  (4.10)) 

Relativistic  particle  momentum;  also  a  transformation 
parameter  in  the  conversion  from  B~L  to  R-\  (Ref.  P3, 

Eq.  (1.13)) 

Counter  symmetry  axis  unit  vector 

Component  of  particle  momentum  perpendicular  to  the 
magnetic  flux  density  vector  B 

Component  of  particle  momentum  parallel  to  the  magnetic 
flux  density  vector  B 

Relative  spectral  intensity  of  frequency  f  (Ref.  PI, 

Eq.  (Al.l)) 

Distance  of  a  particle  from  the  center  of  the  earth  (km) 

Vector  from  center  of  earth  to  tracking  station  in 
Vernal  Equinox  coordinates  (km) 


Pl-5 


Vector  from  center  of  earth  to  satellite  in  Vernal 
l.quinox  coordinates  (Km) 

.  7 

Radius  of  the  •’arth,  0.63712  x  10  meters 

Equatorial  distance  from  dipole  center  to  B  line  (km) 

Extrapolated  range  of  particles  penetrating  an  absorber 
(Ref.  PI,  P3,  P4;  gm  cm~2) 

Range  for  50  percen^  transmission  of  radiation  through 
an  absorber  (gm  cm'  ) 

Scaling  factor  -  twice  the  number  of  accumulated  particle 
counts  corresponding  to  one  frequency  change  on  the  data 
record  (particle  counts  per  flip-flop) 

Normalized  particle  momentum  (Ref.  P3,  Sec.  4.23) 

Detector  A  response  function  (Ref.  P4) 

Matrix  element  representing  the  kernel  of  the  integral 
equation  relating  the  true  count-rate  C(0)  and  the 
directional  flux  density  -1(0') 

A  parameter  in  the  transformation  from  B-L  to  R-X 
(Ref.  P3,  Eq.  (1.14)) 

Time  interval  selected  for  digitazation  (usually  24  sec) 
Longitudinal  drift  period  (Ref.  P3,  Eq.  (3.2S)) 

Bounce  period  (Ref.  P3,  Sec.  3.26)  (sec) 

Cyclotron  period  (Ref.  P3,  Eq.  (3.9))  (sec) 

Bounce  pe-*-4  id  of  a  particle  with  mirror-point  at 

latitude  defined  by  Eq.  (2.8)  (Ref.  P2)  (sec) 

pin 

Time  at  which  particle  count-rate  for  interval  between 
i^iquency  cnange  (flip)  is  assigned  (center  of  interval) 
(UT  min) 

Initial  time  of  time  interval  T  (UT  min) 

Vector  from  tracking  station  to  satellite 
Alfven  velocity  (Ref.  P3,  Eq.  (4.6))  (km  sec-*) 

Speed  of  particle 

Total  energy  (Ref.  P3,  Sec.  3.2) 

Axis  of  maximum  moment  of  inertia  Explorer  4 


Axes  of  Vernal  Equinox  inertial  coordinate  system;  X 
in  direction  of  Ve_rnal  Equinox  of  1951;  Y  in  earth's 
equatorial  plane;  Z  to  form  right  hand  orthogonal  system 

Axes  of  satellite  fixed  rotating  coordinate  system;  X 
space  fixed  in  direction  of  angular  momentum  vector  L; 

Yy  in  the  direction  of  satellite  symmetry  vector  A 
(rotating  in  a_j>lane  perpendicular  to  L  in  propeller¬ 
like  motion) ;  Zf  to  form  a  right-hand  orthogonal  system 

Axes  of  satellite  inertial  coordinate  system;  X  in 
direction  of  angular  momentum  vector  L;  X  -Y  plane 
oriented  to  include  B  the  magnetic  flux  lenlity  unit  vector 
at  time  zero  when  an  r-f  null  is  observed;  Z  to  form  a 
right-hand  orthogonal  system 

First  order  Bessel  function  irregular  at  the  origin 

Atomic  number  of  scattering  atoms 

Axis  of  minimum  moment  of  inertia  Expioi-er  4 

Pitch  angle.  The  angle  between  the  particle  momentum 
vector  and  the  magnetic  flux  density  vector  (deg) 

Pitch  angle  of  the  particle  that  defines  the  loss  cone 
(deg) 

Ratio  of  particle  speed  to  speed  of  light  in  free  space 

Relativistic  mass  ratio  (Ref.  P3), 

Spectral  exponent  (Ref.  P3,  Eq.  (4.19)) 

Phase  angle  of  tumble  motion  of  satellite 

Geodetic  latitude  (radian) 

Geocentric  latitude  (radian) 

Geocentric  latitude  of  satellite  (radian) 


Phase  angle  of  roll  (spin)  motion  of  satellite 


Detection  efficiency  of  counter 


r  d, 


( AQa ^ 


Hour  angle  of  the  Vernal  Equinox  (radian) 

Angle  between  the  Greenwich  meridian  and  Vernal  Equinox  at  t 
East  longitude  of  the  tracking  station  (radian) 


t  Fuler  an^le  specifying  the  opening  angle  of  precession 

about  angular  momentum  vector  L  (Kef.  PI,  P4,  Fig.  2,b) 
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e* 
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^•gm' 
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t.  T 

i>  V 
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ifr(e') 
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coP 


Angle  between  the  directional  counter  symmetry  axis  and 
the  plane  perpendicular  to  the  magnetic  flux  density  B; 

0  =  n/2  -a  (deg);  A  =  n/2  -0 

Angle  between  directional  flux  density  vector  J  and  the 
plane  perpendicular  to  the  magnetic  flux  density  vector 
B  (deg) 

East  longitude  of  satellite  fr^m  Vernal  Equinox  (radian) 
Geographic  latitude  (radian) 

Conjugate  geographic  latitude  (radian) 

Geomagnetic  latitude  (radian) 

Magnetic  dipole  moment 
Time  parameter 

G-M  counter  deionization  (deadtime)  (sec) 

Scintillation  counter  deadtime 
Whistler  frequency  (Ref.  P3,  Eq.  (4.13)) 

Cyclotron  frequency  (Ref.  P3,  Eq.  (4.2)) 

Hm  wave  frequency  (Ref.  P3,  Eq.  (4.2)) 

Flux  invariant  (Ref.  P3,  Eq.  (3.28)) 

Euler  angle  measured  in  X-Y  plane 

Pi'ecession  rate  about  angular  momentum  vector  L  (rad  sec'*) 

Euler  angle  in  the  plane  of  rotation  about  Z'-axis  from 
the  line  of  nodes  to  the  rotating  satellite  fixed  X'-axis 
(deg)  . 

Element  of  an  orthonormal  set 

Residual  angular  velocity  about  Z'-axis  after  transition 
to  propeller- like  configuration  (rad  sec'  ) 

Plasma  frequency  (Ref.  P3,  Eq.  (3.28)) 

Doppler  shifted  whistler  frequency  (Ref.  P3,  Eq.  (4.7)) 

Angular  velocity  about  Z'-axis  (rad  sec"*) 

Initial  angular  velocity  about  the  Z'-axis  at  orbit  in¬ 
sertion  (rad  sec  ) 


Pi-6 
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A.  Summary  of  Pn.D.  Dissertation,  Donald  J.  Manson,  "Van  Allen 

Radiation  Belt  and  ARGUS  Directional  Flux  Density  Distributions, 
Explorer  4,  Satellite  Data1';  Saint  Louis  University,  Physics 
Department  (February  1967). 

Al.  Introduction 

A  recent  determination  (Ref.  Al)  of  the  correct  roll  (or  spin) 
rate  for  satellite  Explorer  4  has  allowed,  for  the  first  time,  system¬ 
atic  and  reliable  determinations  of  directional  flux  density  distribu¬ 
tions  of  geomagnetically  trapped  charged  particles  Injected  by  the 
ARGUS  Event  1  of  27  August  1958.  The  results  of  the  present  analysis 
(1)  indicate  a  pronounced  and  unique  omnidirectional  penetrating  radia¬ 
tion  effect  associated  only  with  ARGUS  Event  1,  (2)  include  a  comparison 
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of  directional  flux  density  distributions  (particles  cm  sec  ster  ) 
of  the  naturally  (inner  V-n  Allen  Belt)  and  artificially  (ARGUS  Event  1) 
injected  radiation,  and  (5)  permit  determination  of  a  decay  curve  for  a 
limited  time  period  of  the  measured  maximum  (perpendicular  to  B)  direc¬ 
tional  particle  flux  density  of  the  artificially  injected  (ARGUS  Event  1) 
radiation. 

A2.  Analysis  and  Experimental  Details 

The  angle  6  between  the  directional  (scintillation  Detector  A, 
Channel  2,5)  counter  arid  a  plane  perpendicular  to  the  magnetic  flux 
density  vector  B  is  determined  as  a  function  of  time  by  calculating 
the  rotational  (roll  and  tumble  motion)  and  forward  motions  of  the 
satellite  i.»  the  Vernal  Equinox  coordinate  system.  Critical  phase 
angles  for  the  roll  and  tumble  motions  at  an  initial  reference  time 


are  determined  by  analysis  of  the  r-f  electric  field  strength  patterns 
from  the  telemetry  and  by  "loop  closure-search"  (based  on  the  assumption 
that  at  equal  angles  of  the  counter  axis  relative  to  B  the  counting  rates 
are  equal;  see  Ref.  A2  for  details).  Each  short  time-period  section  of 
count-rate  data  between  signal  fadeouts  (corresponding  to  one-half  a 
tumble  cycle)  is  analysed  separately.  Satellite  telemetry  yields  pri¬ 
marily  count-rates  versus  time  (Ref.  A2)  which,  combined  with  the  cal¬ 
culation  of  0  as  a  function  of  time,  yields  observed  count-rate  as 

a  function  of  0.  Correcting  for  deadtime,  counter  efficiency  and  effec¬ 
tive  area  (Ref.  A2),  and  solving  the  integral  equation  (Ref.  A3)  relating 
count-rate  to  the  directional  flux  and  the  counter  response  function, 

yields  by  an  unfolding  procedure  the  directional  flux  density  J(0‘) 
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(particles  cm  sec  ster  ).  With  the  aid  of  the  magnetic  mirror  con¬ 
dition  J  „(0'),  (0',  angle  between  vector  J  and  a  plane  perpendicular 

ulclX 

to  B),  the  directional  particle  flux  density  perpendicular  to  B  is  then 
calculated  for  larger  values  of  B  of  the  L  shell  on  which  the  J(0')  vs 
0'  distribution  data  was  acquired  by  Explorer  4. 

Figure  A1  demonstrates  the  validity  of  the  method  used  to  calcu¬ 
late  0  as  a  function  of  time  by  indicating  the  consistency  of  all  the 
data  acquired  by  the  directional  counter  for  4  minutes  of  trajectory 
in  the  Natural  region  of  radiation  (L  =  1.35-1.44;  approximately  1200 
miles  of  satellite  trajectory).  It  is  evident  that  much  of  the  observed 
scatter  is  due  to  decrease  in  the  maximum  count-rate  as  the  satellite 
moves  to  increased  values  of  B  and  L,  and  that  the  shape  of  the  count- 
rate  vs  0  distribution  is  sensibly  constant  with  small  variations  of  L 
for  the  Natural  radiation  in  this  region.  The  relatively  constant 
counting  rate  for  large  values  of  0  indicates  a  penetrating  component 
of  the  radiation. 


A2 


A3.  Results 


All  the  data  acquired  from  the  unshielded  G-M  counter  (Detector  C, 
Channel  3)  and  the  directional  scintillator  (Detector  A,  Channel  2)  for 
Passes  4l4  and  427  are  summarized  in  Figs.  A2  (Upper  and  Lower).  Tri- 
angle-points  for  each  of  the  two  Figures  are  count-rate  maxima  (counter 
axis  perpendicular  to  ll)  determined  from  directional  distributions  such 
as  shown  by  Fig.  Al,  dot-points  are  in  turn  the  minima  (counter  axis 
parallel  to  1') ;  the  dashes  are  the  average  count-rate  from  the  unshielded 
G-M  counter  (detects  electrons  of  energy  greater  than  3  MeV;  protons, 
greater  than  31  MeV)  and  the  X-points  are  the  average  count-rates  from 
the  shielded  (Detector  D,  Channel  1)  omnidirectional  G-M  counter  (detects 
electrons  of  energy  greater  than  6  MeV;  protons,  greater  than  40  MeV). 
Decrease  in  maximum  count-rate  with  increasing  B  and  L  (Fig.  A2)  is 
evident  for  the  Natural  radiation  (to  the  "left"  of  the  arrows  along 
the  time  axes);  clearly  visible  is  ARGUS  Event  1  at  L  =  1.62  -  L  =  1.7 6. 
For  Pass  4l4  a  core  of  penetrating  radiation  is  indicated  (between  arrows 
along  time  axes)  at  L  =  1.71  to  1.76,  1.8  hrs  following  the  nuclear  event; 
this  directional  penetrating  component  is  .ot  observed  24  hrs  later  {Pass 
427).  Comparison  of  Passes  414  and  427  indicates  essential  constancy 
of  the  Natural  radiation  but  rapid  decay  of  the  ARGUS  shell  within  a 
24-hour  period,  and  also  indicates  a  change  in  energy  distribution 
suggested  by  a  difference  in  ratios  (of  the  penetrating  component  of 
the  directional  scintillator  count-rates  to  the  average  count-rates  of 
the  omnidirectional  unshielded  G-M  counter)  for  the  same  24-hour  period. 

Figure  A3  indicates  a  comparison  of  the  directional  flux  density 
distribution,  J ( 0 1 )  vs  4\  for  the  Natural  (left-hand  plot),  ~nd  the 
ARGUS  Event  1  radiation  (right-hand  plot)  of  28  August  1958,  25.48  hrs 


following  the  event. 


A4.  Short-lived  Penetrating  Radiation  Phenomena  Possibly  Related  to 
Geomagnet  lea liy  Trapped  Debris 

This  section  contains  a  revised  (relative  to  our  Final  Report, 
DASA  Report  l6l3,  15  January  1965)  interpretation  of  the  penetrating 
radiations  observed  during  Pass  4l4  (over  Huntsville)  of  Explorer  4. 
Hence  Figs.  A4,  5.  6  from  Report  l6l3  are  reproduced  in  the  present 
section. 

Figures  A4,  A5,  and  A6  are  the  result  of  detailed  directional 
analysis  of  the  data  acquired  from  the  region  Indicated  by  the  arrows 
on  time  axis  (Fig.  A2,  Upper)  where  the  penetrating  radiation  mentioned 
in  A3  above  was  observed.  The  method  of  treating  each  limited  portion 
of  telemetry  data  between  signal  fadeouts  allows  plotting  of  the  count- 
rate,  vs  time  and  0,  during  each  fractional  part  of  a  rotation  of  the 

mti |» 

scintillator  counter  axis  with  respect  to  B.  As  indicated  by  Fig.  A1 
this  method  yields  a  consistent  disk-like  shape  for  the  Natural  radia¬ 
tion  with  all  data  clustering  around  one  curve. 

Figures  A4,  A5,  and  A6  indicate  considerable  departure  of  the 
ARGUS  Event  1  radiation,  1.8  hrs  after  the  event,  from  the  disk-like 
shape  of  the  Natural  radiation  and  so  each  partial  rotation  of  the 
counter  axis  with  respect  to"lT  has  been  plotted  separately  except  for 
Fig.  A4  (b,c,d)  where  three  such  rotations  (from  the  early  penetration 
of  the  ARGUS  shell  have  been  superimposed  to  yield  a  count-rate  vs  0 
distribution  much  like  that  for  the  Natural  radiation). 

In  Fig.  A4(f)  the  time  sequence  shown  by  the  arrows  indicates 
Inal  the  ave  -age  count-rate  is  increasing  rapidly  during  one  partial 
rotation  (approx.  3.6  secs)  so  that  at  equal  angles  of  the  counter 
axis  with  respect  to  B  (e.g.,  0  -  oO  deg)  two  different  count-rates 


are  observed,  the  later  count-rate  having  increased  as  the  satellite 
penetrated  farther  into  the  ARGUS  shell.  It  should  be  noted  that 
the  omnidirectional  penetrating  component  (at  large  0  values)  has 
increased  rapidly  (from  2000counts  sec  1  in  Fig.  A4(b,c,d)  to  approxi¬ 
mately  5000  counts  sec"1  in  Fig.  A4(f).)  The  large  count-rate  (9000 
sec"1,  0  =  92  deg,  Fig.  A4(f)),  at  larger  angles  of  0  indicates  a 
departure  from  the  disk-like  shape  of  the  radiation.  Figure  A4(g) 
indicates  that  the  scintillation  counter  is  saturating  (count-rate, 
greater  than  10,000  sec"1)  when  perpendicular  to]?  (0  =  0)  while  the 
count-rate  at  0  =  50  deg  has  dropped  (from  9000  to  7000  sec  1). 

Fir-re  A4(h)  indicates  saturation  at  small  to  moderate  0  angles 
(  >  10,000  sec"1,  0  40  deg)  and  a  count-rate  of  almost  7000  sec  1 

at  0  =  60  degrees  -  again  a  considerable  departure  from  disk- like 
shape  of  the  count-rate  vs  0  distribution. 

Figure  A5(j)  (acquired  2.5  secs  later  than  Fig.  A4(h))  indicates 
a  more  normal  disk-like  shape  but  for  a  limited  range  around  0=0. 
Figures  A5(j,k,l  and  ra)  if  superimposed  would  yield  a  single  consistent 
distribution  much  like  Fig.  A5(m)  where  the  Cobg  vs  0  distribution  is 
disk-like  (0  £  30  deg)  but  with  a  large  omnidirectional  penetrating 
background  component  (  ^  5000  sec  1).  The  final  5  points  of  the  time 
sequence  (arrows)  shown  in  Fig.  A5(m)  indicate  a  sudden  saturation  of 
the  counter  (0  £  30  deg). 

This  saturation  is  continued  in  Fig.  A6(n)  at  0  =  30-35  degrees 
but  is  followed  in  Figs.  a6(o)  by  an  abrupt  change  to  a  more  nearly 
disk-like  distribution;  followed  immediately  in  turn  by  saturation  for 
all  angles  0  ^  40  deg.  Figure  a6 (p )  (0.2  sec  later)  indicates  a 

return  to  the  normal  disk-like  shape  followed  by  a  sharp  drop  in  the 


count-rate  at  6  -  0  as  the  satellite  leaves  U.a  ARGUS  shell.  Figure 
At (q)  is  for  the  last  measurable  count-rates  from  the  ARGUS  Event  1 
wing  as  the  satellite  moves  into  the  slot  between  the  two  Van  Allen 
belts  away  from  the  ARGUS  shell. 

The  analysis  presented  relative  to  Figs.  A4,  5.  6  suggests  that 
the  data  from  Pass  414,  1.8  hrs  after  Event  1  contains  evidence  for 
short  term  penetrating  components  (directional  and  omnidirectional) 
of  the  radiation  which  are  neither  spatially  homogenous  nor  Isotropic. 
The  fact  that  the  scintillation  counter  goes  to  saturation  suddenly 
during  one  portion  of  the  rotation  of  the  counter  axis  with  respect 
to  vector  B  as  in  Fig.  A5(m)  and  A6(o)  while  the  remaining  portion 
of  the  rotation  measures  a  fairly  normal  distribution  suggests  a  vary¬ 
ing  directional  component  of  the  flux  at  these  points  (possibly  due  to 
gamma  radiation  from  trapped  debris  near  the  ARGUS  shell).  The  normal 
disk-like  distributions  (Fig.  A3)  measured  in  the  ARGUS  shell  1  day 
later  for  Pass  427  apparently  are  present  in  Pass  4l4  but  are  obscured 
and  distorted  by  the  large  penetrating  omnidirectional  component  and 
bursts  of  directional  radiation  which  saturate  the  counter.  These 
effects  are  greatly  diminished  in  Pass  415,  (3.67  hrs  following  Event  1) 
and  are  not  found  in  Pass  427  (25.48  hrs  following  Event  1). 

No  such  effects  were  found  from  ARGUS  Event  2  where  data  from 
Pass  454  (I.87  hrs  following  Event  2,  30  August  1938;  Refs.  A2,  A4) 
indicate  the  usual  disk-like  distributions. 

Figure  A7  indicates  the  time  decay  of  maximum  directional  flux 
density  f°r  two  d&ys  following  ARGUS  Event  1.  It  was  not 

possible  to  follow  the  decay  for  a  longer  period  than  »,wo  days  because 
by  the  third  day  tne  directional  trapped  radiation  intensity  has  dimin¬ 
ished  to  near  background  value. 

Ar, 


The  x-points  lie  well  on  a  straight  line  (exponential  decay  law), 
the  solid-line  curve  is  a  plot  of  the  t  law  (observed  for  the 

omnidirectional  radiation,  see  Part  C  of  the  present  report)  for  com¬ 
parison.  Because  of  the  snort  time  Interval  represented,  about  2 
days,  a  decision  between  the  two  decay  laws  is  not  considered  valid. 

Figure  A8  is  a  3-dimensional  plot  of  J  (0'),  the  directional 
°  max 

flux  density  perpendicular  to  vs  B  and  L;  both  Natural  and  ARGUS 
Events  1  and  2  radiations  are  shown.  Figure  A8  summarizes  most  of 
the  data  analysed  for  the  scintillation  counter,  Detector  A. 

A5.  Conclusions 

Detailed  analysis  of  ARGUS  directional  flux  density  distributions 
from  Explorer  4  data  (Passes  414  and  427,  over  Huntsville,  Alabama,  on 
27  and  28  August  1958  respectively)  for  1.8  and  25.48  hrs  respectively 
after  Event  1  and  comparison  with  Natural  flux  density  distributions 
indicate  that  the  disk-like  shape  as  measured  by  the  scintillation 
counter  in  the  ARGUS  shell  is  not  substantially  different  from  the 
Natural  radiation  (by  Fig.  A3,  the  half -width  at  half-maximum  fhwhm) 
is  7  deg  broader  on  the  average  for  the  ARGUS  J(0')  compared  to  the 
Natural  radiation;  maximum  values  of  directional  flux  densities  are 
12  x  107  for  the  Natural  vs  6.7* 107  cm  ^sec  ^ster"^  for  ARGUS).  The 
27  August  1958  data  from  Pass  4l4  (Ref.  A2)  and  Pass  415  (Ref.  A3) 
indicate  that  these  measured  distributions  are  temporarily  obscured 
and  distorted  by  a  large  penetrating  omnidirectional  component  and 
a  varying  directional  component.  The  short-time  limited  decay  data 
suggest  an  exponential  decay  of  maximum  directional  flux  density  at 
L  =  1.7  but  a  t  law  observed  for  the  omnidirectional  data  from 


G-M  counter  (Detector  C)  Is  not  excludod.  It  ruf.y  be  suggested 


that  particle  and  energy  differences  between  Natural  and  ARGUS  radia¬ 
tions  can  be  expected  to  yield  somewhat  different  decay  laws. 
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Pig.  Al.  Pass  427  (over  Huntsville,  Alabama,  28  August  1958) 

observed  directional  count -rate  C  .  vs  0  for  4  min 

obs 

in  the  region  of  Natural  radiation. 

o,  Cobg  vs  0;  227-5  -  229-5  min  UT;  L,  1.35  -  1-42 
earth  radii;  B,  0.14^  -  0.164  gauss;  ,  CQbs  vs  0; 
229.5  -  231.5  min  UT;  L,  1.42  -  1.49  earth  radii; 

B,  0.164  -  .180  gauss. 

Data  were  plotted  by  combining  about  one-half  of 
the  100  individual  sections  of  CQbs  vs  0  described  in 
Ref.  A2,  Sec.  3.3.  All  data  acquired  within  the  time 
period  227.5  -  231.5  UT  are  plotted  to  show  the  amount 
of  scatter  in  the  data  and  to  demonstrate  the  consist¬ 
ency  of  the  analytic  interpretation  using  the  experi¬ 
mentally  determined  parameters:  roll  and  tumble  phase 
angles  8  and  V  ,  and  the  roll  and  tumble  rates. 
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o 


PASS 427  NATURAL 
227.5-231.5  UT(min) 
AUGUST  28,1958 


THETA  (deg) 


Pig.  A2.  (Upper,  Lower) 

Natural  and  ARGUS  Event  1  count-rates  vs  time  and 
B-L  coordinates  for  Passes  4l4  and  427  over  Huntsville, 
Alabama.  Comparison  of  Average  Omnidirectional  Count- 
Rate  with  Maximum  (perpendicular  to’B)  and  Minimum 
(parallel  to"£)  Directional  Counting  Rates. 

Triangle-points,  maxima;  and  dot -points,  minima, 
scintillation  counter  (Derector  A)  observed  count-rates; 

Dash-points,  15-sec  averages  of  observed  omnidirectional 
(unshielded  G-M  counter.  Detector  C)  count-rates; 

X-points,  6-sec  averages  of  observed  omnidirectional 
count-rates  (shielded  G-M  counter,  Detector  D). 

Arrows  on  the  time  axis  indicate  time  region  in  which 
distributions  of  Pigs,  A4,  5.  and  6  (Pass  4l4)  were  acquired 
(ARGUS  Event  1).  Calibration  information  for  the  three 
counters  is  summarized  in  the  Preface  (PI)  of  the  present 
report . 
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Pig.  A3.  Directional  flux  density  distributions  J(0')  vs  0 
for  Pass  427,  28  August  1958;  Natural  (left-hand 
plot)  and  ARGUS  Event  1  (right-hand  plot)  radiations. 
Variations  of  magnitudes  of  *Jmax(0‘)  =  J(0)  and  shape 
of  the  distributions  are  shown.  J(0')  for  ARGUS 
Event  1  radiation  were  corrected  for  varying  penetrating 
radiation  component  as  explained  in  Ref.  A2. 


All 


Pigs.  A4,  A5,  and  A6 

Detailed  sequence  of  scintillation  counter  count-rates 
Cq^s  vs  ®  (angle  between  scintillation  counter  axis 
and  the  plane  perpendicular  to  *B)  for  Pass  414  ARGUS 
Event  1  penetration,  27  August  1958*  Where  used 
arrows  indicate  time  sequence  in  acquisition  of  points. 
Distributions  indicate  a  large  omnidirectional  pene¬ 
trating  component  and  a  varying  directional  component 
of  the  radiation  (discussed  in  detail  in  the  text). 
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ARGUS  Event  1 


ARGUS  Event  1 


ARGUS  Event  1 


Pig.  A7.  Decay  of  maximum  directional  particle  flux  densities 

J  ( 0 * }  vs  time  for  ARGUS  Event  1.  Plotted  points 
max 

are  for  flux  density  distributions  J(0*)  at  0'=  0. 

The  maximum  ARGUS  Event  1  radiation  is  observed  to 
decay  one  order  of  magnitude  in  45  hrs.  The  lower 
curve  is  a  plot  of  the  t  law  (observed  for  the 

omnidirectio- ' 1  radiation,  Part  C  of  the  present 
report,  for  the  same  ARGUS  Event  1  radiation). 
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(crtT^sec"1  ster”* 


DECAY  OF  MAXIMUM 
OMNIDIRECTIONAL  FLUX  DENSITY 


FROM  ARGUS  EVENT  1 


_j _ I _ i _ L  -  i  I  i  I  ■  I  ■  I 

10  20  30  40  50  60 

TIME  (hrs.  after  time  zero,  event  I) 


Pig.  A8.  Directional  maximum  particle  flux  density  Jmax(QI) 
vs  B  and  L  for  both  Natural  and  ARGUS  events  1  and 
2  radiations  (Detector  A,  Channel  2) 


Detailed  information  is  tabulated  as  follows. 


Natural  Radiation 


L  Slope  Pass  No.  Date 

(earth  radii) 


oL 

p 

1.360 

r .  . . . 

-7.92 

28  Aug. 

1938 

1.380 

-4.82 

27  Aug. 

1938 

% 

1.407 

-7.15 

427 

28  Aug. 

1958 

o 

1.409 

-2.76 

414 

27  Aug. 

1958 

£ 

1.434 

-7.13 

454 

30  Aug. 

1958 

& 

1.458 

-7.19 

454 

30  Aug. 

1958 

7 

1.580 

-7.92 

427 

28  Aug. 

1958 

ARGUS  EVENT  1 


e 

1.690 

-3.78 

414 

27  Aug.  1958 

Z 

1.694 

-6.42 

427 

28  Aug.  1958 

H 

1.724 

-4.10 

427 

28  Aug.  1958 

F 

1.728 

-4.11 

440 

29  Aug.  1958 

■? 

1.730 

-3.58 

415 

27  Aug.  1958 

$ 

1.743 

-5.05 

427 

28  Aug.  1958 

0 

1.748 

-2.00 

427 

28  Aug.  1958 

X 

1.777 

-4.48 

427 

28  Aug.  1958 

ARGUS  EVENT  2 


IT 

2.124 

-7.55 

454 

30  Aug.  1958 

f> 

2.124 

-3.60 

454 

30  Aug.  1958 

cr1 

2.144 

-15.40 

454 

30  Aug.  1958 

r 

2.172 

-10.40 

454 

30  Aug.  1958 
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B.  Summary  of  M.3.  Thesis  of  Joseph  F.  Fennell,  "Van  Allen  Belt  and 


ARGUS  Directional  Flux  Density  Distributions,  Explorer  4  Satellite 
Data,  ARGUS  Event  2",  Saint  Louis  University ,  Physics  Department 
(July  1966). 

Bl.  Introduction 

This  Thesis,  Reference  Bl,  contains  essentially  an  analysis  of  the 

-2  -1  -1 

directional  distribution  of  flux  densities  (particles  cm  sec  ster  ) 
vs  the  angle  between  the  counter  symmetry  axis  and  the  plane  perpendi¬ 
cular  to  the  magnetic  flux  density  vector,  of  geomagnetically  trapped 
naturally  occurring  and  artificially  injected  charged  particles  asso¬ 
ciated  with  ARGUS  Event  2  for  Pass  454  over  Huntsville,  Alabama  (for  a 
time  interval  of  about  10  min.  beginning  1.8  hrs  following  the  Event). 

The  purpose  of  this  limited  analysis  was  to  compare  the  situation  with 
that  of  Pass  4l4  (Manson  Dissertation,  Ref.  B2)  with  the  Passes  occur¬ 
ring  approximately  3  days  apart.  The  count-rate  data  were  obtained 
from  a  scintillation  photomultiplier  detector  (Detector  A,  Channel  2 
Ref.  B2) ;  window  aperture,  15  degrees  half-width  at  half-maximum; 

Agff£,  0026  cm  ;  deadtime  100  microsecs  aboard  Satellite  Explorer  4. 

Tne  thesis  represents  the  first  detailed  examination  of  the  direc¬ 
tional  data  for  ARGUS  Event  2  by  the  Saint  Louis  University  group.  There 
is  presented  a  short  discussion  of  a  method  for  determining  the  proper 
roll  phase  angle  used  in  this  analysis  and  the  results  from  determinations 
of  satellite  roll  frequencies  for  times  at  which  they  had  not  been  pre¬ 
viously  determined  with  sufficient  preclsiu.,  to  yield  reliable  determi¬ 
nations  of  observed  count-rate  vs  an6^e  6  (between  the  counter 

symmetry  axis  and  a  plane  perpendicular  to  the  magnetic  flux  density 
vector  B) . 

Bl 


B2.  Determination  of  the  Phase  Angle  for  the  Roll  Motion  of  Satellite 


To  determine  the  directional  flux  density  J(@')  vs  angle  0'  between 

the  counter  axis  and  the  incident  flux  direction,  C  .  (0)  as  a  function 

obs 

of  0  must  be  known.  9  at  time  t  (relative  to  an  arbitrary  reference  time) 
is  given  by 

6  -  A rc  s  i  n[  b,  S  i nOPt  *8)  -  bxCos  Cpt  +Sj  S  incfit  +Y ) 


+  b3Cos  (^ir  +3)  CosC^t+Y)], 


(Bl) 


where  all  parameters  are  known  explicitly  except  3  and  7^  •  the 

roll  and  tumble  motion  phase  angles.  y  was  determined  as 

7  -  Arccos  (A  •  Ys) ,  <m> 

\ 

where  A  is  the  satellite  symmetry  axis  vector  and  Yg  is  a  unit  vector 

in  Vernal  Equinox  Coordinates  (Ref.  B2). 

3  is  determined  (Ref.  Bl)  in  Eq.  (Bl)  by  assuming  counting 

rates  are  equal  approximately  at  equal  angles  0  in  the  Natural  radiation 

region  (over  short  periods  of  time,  t  <  0.5  min).  Taking  two  times 

t  ,  t,  with  equal  counting  rates,  it  can  be  shown  (Ref.  Bl)  from  Eq.  (Bl) 
*  J 

that 

Sjj  -  Arcta.n  (~fjj  /9jj)  t»3) 


where  f, .  =  f,  +  f . ,  g, ,  =  g,  -*•  6. > 

■» -J  -  u  -J  t  1  J  . 

f  =  b  (t  )  SinClj/tj)  +  P(t  1,y  )  Cos(  xp  ti), 

*  • 

£i  =  bjUj)  cosn/ztj  -  PUj.y )  sm{  xj;  ti) 

p(tA ,  y )  =  b  (t1)  Cos ( <£>  t1+  y  )-b2(t1)  Sin(  <$>  t^y). 


and 
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An  average  <j  is  then  determined  by  considering 

§a.v<]  =  A  ret  a_n  ("  L  (t/j  /9ij ) ) 

N 


(B4) 


where  the  sum  is  over  N  time  pairs  t^,  tj.  The  consistency  of  the  data 
can  be  seen  by  considering  the  variations  in  the  §  .y  from  5  avg. 

If  the  spread  in  the  values  of  S  . .  is  small  the  random  fluctuations 

s 

in  the  data  are  small.  ^ivg  was  determined  (with  respect  to  an  arbi¬ 
trary  reference  time)  as  94°  for  Pass  454  and  this  value  reproduced  the 
count-rate  variations  in  a  simulator  program  (see  Ref.  B2). 

The  roll  rate  for  Explorer  4  was  determined  for  the  entire  tele¬ 
metry  times  of  Passes  4l4  through  454  by  fitting  a  least  square  poly¬ 
nomial  to  the  known  roll  rate  values  at  discrete  reference  times  of 
Passes  414,  415,  427,  and  454.  A  second  method  for  determining  the  roll 

rate,  based  on  the  fact  that  Sin  0  is  a  minimum  when  C  .  (0)  is  maximum 

obs 

is  presented  by  Fennell  (Ref.  Bl).  (The  second  method  was  applied  to 
Pass  440  and  the  resulting  roll  rate  was  used  by  J.  M.  Paikeday  (Ref.  B3) 
to  obtain  a  directional  flux  density  distribution  for  Pass  440  ARGUS 
Event  1  radiation. ) 


B3.  Results 

With  0  determined  as  a  function  of  time,  plots  CQbs(0)  vs  0  are 
made  and  then  corrected  for  the  counter  deadtime  and  (effective  area 
x  counter  efficiency)  to  yield  the  true  count-rate  C(0) 


C(9)  = 


Co  bS  (&) 


(B5) 


eAfiff  (1  rT>  Cobs  (6)  ) 


The  corrected  or  true  count-rates  0(6)  are  used  to  compute  the  corres¬ 
ponding  directional  flux  densities  J  ( 6 1 )  as  described  by  J .  M.  Paikeday 
(Ref.  B3).  Representative  examples  of  the  observed  count-rates  CQ^s(6) 
vs  0,  true  count-rates  C(0)  vs  0,  and  directional  flux  densities  J ( 0 1 ) 
vs  0',  for  Pass  454  Natural  radiation  are  shown  by  Pigs.  B1  and  B2.  The 
CobS(0)  vs  0  curves  showed  an  "open  loop"  phenomena  (see  Pig.  B3(a)) 
which  is  due  in  part  to  the  rapid  rate  of  change  of  both  the  omnidirec¬ 
tional  and  directional  components  of  the  radiation  as  the  satellite  moves 
through  the  ARGUS  Event  2  shell.  The  curve  was  reduced  to  a  single  valued 
curve  by  subtracting  out  this  component  (compare  Pig.  B3(a)  and 
of  Pig.  B3(c)).  The  resulting  CQ^s(0)  vs  0,  C(0)  and  J(6')  vs  0,  0' 
curves  for  the  ARGUS  radiation  are  shown  in  Figs.  33,4.  Figures  Bl, 2,3,4 
and  Tables  Bl,  B2  in  the  Pigs.  B2,  B4  captions  represent  the  results  of 
this  analysis. 

B4.  Conclusions 

The  maximum  directional  flux  density  and  observed  count-rate  vary 
considerably  in  the  six  minutes  of  the  Natural  radiation  observed  for 
Pass  454  while  the  half -width  at  half-maximum  (hwhra)  remains  relatively 
constant  at  17  +  5  degrees  (see  Table  Bl,  Pig.  B2  caption).  The  varia¬ 
tion  in  maximum  count -rate  and  flux  density  is  probably  related  to  the 
large  spatial  and  magnetic  flux  density  changes  involved  in  the  satellite 
motion. 

The  non-zero  J(0')  for  0'  50  bo  60  degrees  (Pig.  B2),  in  the 

Natural  radiation  region  is  due  probably  to  pene  rating  (independent  of 
01 )  radiation  capable  of  transmission  by  the  counter  shielding  IE.  3  Mev) ; 
it  is  the  background  count-rate  for  the  detector  in  the  region  of  Natural 
radiation. 
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The  flux  densities  and  hwhm  vary  considerably  for  the  ARGUS  Event  2 
Pass  454  radiation.  By  Pig.  B4  the  ARGUS  radiation  appears  generally  to 
have  a  smaller  hwhm  (16°  vs  l8°)  and  a  larger  maximum  flux  density 
(6.4  x  107  vs  4.9  x  107  particles  cm  ^sec-1ster"*)  than  the  Pass  4^4 
Natural  radiation.  It  is  noted  that  there  are  large  differences  in  B 
and  L  coordinates  between  the  Natural  and  ARGUS  radiations.  The  results 
for  Pass  454  are  included  in  Table  B2  for  comparison  with  previous  analy¬ 
sis  (Refs.  B2,  B5,  B4,  35). 

The  maximum  flux  intensity  and  hwhm  of  the  Pass  454  ARGUS  Event  2 
are  considerably  lower  (6.4  x  10'  vs  2.3  x  10°  particles  cni  sec^ster-1) 
and  narrower  (16°  vs  2lf°)  respectively  than  those  of  Pasa  4l4  for  ARGUS 
Event  1  with  both  Passes  occurring  about  1.8  hours  after  their  respective 
events.  Assuming  both  nuclear  devices  were  the  same,  the  differences  in 
Pass  414  and  Pass  454  ARGUS  results  may  be  due  to  spatial  or  geomagnetic 
affects.  The  data  occurred  at  the  two  L's,  I.72  (Event  1)  and  2.1 
( Event  2). 

Some  tentative  conclusions  can  be  offered  for  Pass  454  Event  2,  as 
follows : 

(1)  ARGUS  radiation  is  more  intense  than  the  Natural  as  is  evidenced  from 
the  Tables  B1  and  B2  (see  Pigs.  Bl,  B4  captions). 

(2)  ARGUS  radiation  appears  to  be  slightly  more  disk-like  than  Natural 
as  is  evidenced  from  the  hwhm  (Tables  Bl,  B2,  B3). 

.3)  ARGUS  radiation  showed  no  sign  of  any  rapidly  varying  omnidirectional 
penetrating  radiation  effect  as  was  observed  by  Manson  (Ref.  B2)  in 
Pass  4)4 
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Since  Pass  4^4  is  the  only  ARGUS  Event  2  directional  data  which 
has  been  analyzed  by  the  methods  presented  In  Refs.  B1  and  B2  to  date, 
a  more  intensive  study  of  Event  2  data  and  a  more  complete  comparison 
with  Event  1  directional  data  must  await  the  analysis  of  more  Event  2 
directional  data. 
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Pig.  Bl.  Representative  °ol3S(®)  vs  ®  curves  for  Pass  4^4  Natural 

radiation.  Actual  data  points  are  shown  to  indicate  the 

scatter  in  C  ,  (0)  vs  0. 
obs 
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Fig.  B2.  True  count-rate  C(0)  (counts  cm  sec  corrected  for  deadtime, 

effective  area  and  counter  efficiency)  and  directional  flux 

-2  -1  -1 

J(0')  (particles  cm  sec  ster  )  vs  0  and  0'  for  the  Fig.  B1 
Cobs(e)  curves.  The  bars  on  the  J(0')  vs  0*  curves  indicate 
the  scatter  in  calculated  values  of  J ( 0 * )  due  in  turn  to  the 
scatter  of  points  in  the  corresponding  CQbs(0)  vs  ®  curves. 


Table  Bl.  Characteristics  of  Natural  radiation  (from  Ref.  Bl) 


Time  (min) 

♦hwhra  (deg) 

J(6') 

max 

B  (gauss) 

L(  earth  radii 

300.00 

22.0 

4.9x10^ 

.1734 

1.376 

301.221 

17.5 

4. 9x10^ 

,1782 

1.414 

501.913 

17.5 

4.7x10; 

.1806 

1.434 

302.370 

17.5 

3.5xlo; 

.1834 

1.458 

304.570 

16.0 

1.4x105 

.1912 

1.583 

305.725 

16.0 

8.2x10^ 

.2172 

1.640 

306.990 

12.0 

7-7x10 

.2316 

1.776 

*  half-width  at  half -maximum 
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C(0)  TRUE  COUNT  RATE  (cm“  sec**  )  jj  C(8 )  TRUE  COUNT  RATE  (cm"2  sec" 
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Pig.  B3. 

(a)  CQbs(e)  vs  0  ARGUS  Event  2  'loop'  before  the  linearly 
increasing  background  has  been  subtracted.  Compare 
Pig.  B3(c)  which  is  the  resulting  CQbs(0)  vs  9  curve 
after  subtraction  of  background.  The  background  sub¬ 
tracted  is  represented  by  (constant  x  t^),  where  t^  is 
the  time  (following  the  time  of  the  first  point  of  a 
given  Cobg(6)  vs  ®  curve)  corresponding  to  the  i°n 
point  of  the  given  curve. 

(b)(c)(d)  c0bS(®)  vs  e  f°r  Pass  ^5^  ARGUS  Event  2  radiation 

after  loop  closure  is  accomplished. 
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20  40  60  20  40  60 
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20  40  60  20  40  60 

0  (deg)  9  (deg) 


Pig.  B4.  C(0)  and  J(0')  vs  0  and  0’  for  the  CQbs(e)  vs  0  curves 
(Pig.  B3),  Pass  454  ARGUS  Event  2.  The  bars  on  curves 
indicate  the  scatter  in  calculated  values  of  J(0')  as 
above  (Fig.  B2). 

Table  B2.  Characteristics  of  the  Pass  454  ARGUS  Event  2  radiation 
(taken  from  Ref.  Bl). 

Time  (min)  hwhm  (deg)  J(0')  v  B  (gauss)  L(earth  radii) 

filaX 


310.178 

12.0 

3.2xl07 

.2720 

2.120 

310.223 

24.0 

5.63x10? 

.2726 

2.124 

310.279 

13.0 

5.93x10? 

.2740 

2.128 

310.336 

14.0 

6.45x10? 

.2742 

2.144 

310.496 

19.0 

2.82X106 

.2764 

2.172 
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C.  Summary  of  Ph.D.  Dissertation,  John  A.  George,  “Omnidirectional 
Fluxes;  Explorer  4  Satellite  Data,  ARGUS  Events  1  and  2”,  Saint 
Louis  University,  Physics  Department  (February  1967). 

Cl.  Introduction 

In  August  and  September  1958,  three  low  yield  nuclear  detonations 
(ARGUS  Events  1,  2,  3)  were  produced  at  high  altitude  over  the  south 
Atlantic  Ocean.  Electrons  emitted  were  subsequently  trapped  in  the 
geomagnetic  field.  The  characteristics  of  this  radiation  were  measured 
by  instruments  on  board  Jason  sounding  rockets  and  in  the  Satellite 
Explorer  4.  A  number  of  papers  have  appeared  in  which  the  results  of 
what  is  known  as  the  ARGUS  Experiment  were  described  (Refs.  Cl,  C2,  C3) . 

The  present  Dissertation  is  concerned  with  further  analysis  of 
data  obtained  by  Explorer  4  for  ARGUS  Events  1  and  2  for  the  G-M 
Detectors  C  (Channel  3)  and  D  (Channel  1);  therefore  data  for  omni¬ 
directional  radiation  is  considered. 

C2.  Experimental  Details 

The  characteristics  of  satellite  Explorer  4  (1958  Epsilon)  have 
been  described  in  detail  by  Van  Allen  (Ref.  C2) .  For  later  calibration 
information,  as  used  in  the  present  report,  see  Preface  P2  of  this 
report.  The  Anton  302  G-M  counters  only  are  involved  in  the  present 
analysis.  A  summary  of  the  characteristics  of  these  two  detectors  is 
given  in  Table  Cl.  These  detectors  are  omnidirectional  and  do  not 
differentiate  the  kinds  of  particles  being  detected.  However  for  the 
ARGUS  radiation,  it  will  be  assumed  that  the  predominant  contribution 
to  the  counting  rates  is  b>  electrons  resulting  from  the  fission 
processes.  Anton  302  G-M  counters  similar  to  those  on  Explorer  4 


Cl 


have  since  been  flown  on  a  number  of  satellites,  and  it  has  been  con¬ 
cluded  (Refs.  C4,  C5,  C6)  that  the  counting  rates  must  be  attributed 
to  penetrating  electrons  of  Ee  >  1.5  MeV.  Motz  and  Carter  (Ref.  C7) 

studied  the  response  of  Anton  302  counters  to  fission-spectrum 

2 

electrons,  They  found  that  up  to  2  gm/cm  of  lead-absorber  thickness 

the  curves  were  as  expected  for  6-transmission;  however,  for  energies 

beyond  this  range  the  transmission  was  enhanced  somewhat  by  bremm- 

strahlung.  Detector  C  (the  "unshielded"  G-M  counter,  Channel  3)  is 

2 

shielded  by  1.2  gm/cm  of  stainless  steel,  which  is  equivalent  to  a  shleld- 

2  2 
ing  of  1.46  gm/cm  of  aluminum.  The  extrapolated  range  Rg  (gm/cm  ) 

versus  energy  E(MeV)  for  aluminum  given  by  Marshall  and  Ward  (Ref.  C8) 

is 


R,.  =  0.528E  -0.094  (Cl) 

b 

For  Detector  C,  Rg  by  Eq.  (1)  corresponds  to  2.95  MeV  electrons 

(Table  C2) ;  for  Detector  D  (the  "shielded"  G-M  counter,  Channel  1) 

2 

which  had  an  additional  shielding  of  1.6  gm/cm  of  lead,  the  total 

2 

shielding  is  equivalent  to  3.89  gm/cm  of  aluminum  and  Rg  corresponds 
to  7.6  MeV  electrons.  For  the  unshielded  Detector  C  at  least  the 
assumption  of  electrons  being  counteJ  directly  seems  quite  valid.  How¬ 
ever,  the  count-rate  of  the  shielded  Detector  D  may  be  enhanced  some¬ 
what  by  bremmstrahlung.  A  plot  of  the  omnidirectional  geometrical 
factors  versus  the  energy  necessary  to  penetrate  the  various  thickness 
of  material,  using  the  information  of  Tables  Cl  an<  C2,  is  given  in 
Fig.  Cl. 

An  automatic  digitizer  was  used  to  analyze  the  telemetry  play¬ 
back  data  and  transform  it  to  punch  cards.  A  computer  program  was 


C2 


written  to  obtain  count-rate  versus  time  for  each  selected  Pass  of  the 
Satellite  fxplorcr  4  over  Huntsville,  Alabama.  The  observed  count- 
rates  CQks  (sec  )  were  corrected  for  the  deadtime  of  the  detectors  and 
associated  electronics  to  obtain  true  count-rates  C  (sec  *)  by  using 
the  approximation  (instead  of  hq.  (Dl) ,  p.  U2  below) 


where  deadtime  is  62.5  microsecs  (Ref.  C2)  and  was  confirmed  in 
the  present  work.  Figure  C2  is  typical  of  the  data  available  for 
analysis. 

The  B  and  L  coordinates  for  the  ARGUS  region  were  obtained  from 
GE  TEMPO  Data  Center  (Santa  Barbara,  California)  through  the  offices 
of  the  Defense  Atomic  Support  Agency.  Figure  C3  shows  the  regions  of 
interest  traversed  by  Explorer  4  in  B-L  space  for  typical  Passes  over 
Huntsville,  Alabama.  Table  C3  lists  the  complete  set  of  Passes 
analyzed  and  are  arranged  in  three  columns,  each  of  which  corresponds 
to  nearly  the  same  region  of  space  relative  to  the  earth.  Tne  orbi¬ 
tal  frequency  of  Explorer  4  was  about  13  orbits  per  day. 

In  order  to  investigate  the  possible  correlation  of  Explorer  4 


ARGUS  count-rates  with  solar  activity. 


3-hr  K 

P 


indices  are  presented 


by  Fig.  C4  for  the  period  under  investigation  27  August  to  4  September 


19S8.  Subsequent  to  the  large  geomagnetic  fluctuations  observed  at 


the  onset,  a  period  of  decreasing  activity  to  relative  magnetic  quiet 


levels  followed.  No  sudden  commencements  were  observed  until  5-5 


September,  a  period  which  is  classified  as  one  of  magnetically  dis¬ 
turbed  days. 


C3.  Results 


ARGUS  Event  1 

The  decay  of  the  fission  produced  electrons  trapped  in  the  earth's 
magnetic  field  is  indicated  in  Fig.  C5,  a  log- log  plot  of  the  true 
count-rate  versus  T  ,  the  elapsed  time  since  Event  1,  at  a  fixed  L 
value.  The  plotted  points  are  the  true  omnidirectional  count-rates 
of  the  unshielded  G-M  counter  (Detector  C,  Channel  3) .  The  count-rates 
are  differences  between  total  and  estimated  background  rates,  and  are 
for  the  414,  427,  ...518  sequence  of  Passes.  The  maximum  intensities 
observed  are  well  represented  by  the  L  =  1.72  shell.  The  peak  location 
did  not  shift,  at  least  within  the  experimental  error,  from  this  L 
shell  during  the  time  interval  indicated.  If  it  is  assumed  C  (sec  *) 
versus  t(hrs)  is  given  by 

C  =  CQt“n,  (C3) 

where  Cq  and  n  are  undetermined  constants,  a  least  squares  fit  of  the 
data  for  L  =  1.72  yields  n  =  1.13,  Cq  =  1.31  x  103;  whence  Eq.  (C3) 
becomes 

C  =  1.31  x  105t'1,13,  0.19  <  B  <0.23.  (C4) 

This  t~*‘*3  law  agrees  quite  well  with  the  decay  rate  determined  by 
Cladis  and  Walt  (Ref.  C3)  from  Jason  rocket  data  for  ARGUS  Event  2 
(L  =  2.1). 

Figure  C6  is  a  plot  of  the  count-rate  ratio  unshielded/shielded 
omnidirectional  G-M  counters.  Data  for  shielded  and  unshielded  G-M 
counter  of  rates  greater  than  10  sec  *  and  30  sec  ^  respectively,  only 
were  used  to  eliminate  data  averaged  over  larger  regions  of  space. 


C4 


The  data  for  Pass  414  may  be  in  considerable  error  due  to  saturation 
of  the  counters.  The  vertical  bars  are  estimates  of  the  spread  of  data 
arising  principally  from  the  scatter  in  the  count-rates  of  the  shielded 
G-M  counter.  The  count-rate  ratios  terminate  at  Pass  466  due  to  the 
difficulty  of  separating  the  ARGUS  electrons  from  the  natural  back- 
ground  radiation.  The  dashed  line  is  the  plot  of  C  =  CQt  *  .  Of 

interest  is  the  sudden  broadening  of  the  ARGUS  Event  1  shell  subsequent 
Pass  453,  Fig.  C7.  The  shell  width  at  half-maximum  is  given  in  terms 
of  ^  L.Pass  453  was  the  first  pass  *0  record  the  effect  of  Event  2 
but  all  data  of  Fig.  C7  are  for  Event  1.  However,  Event  2  occurred 
subsequent  to  the  satellite  traversing  the  Event  1  Pass  453  shell. 

Since  Passes  up  to  492  occurred  during  a  period  of  low  geomagnetic 
activity  (K^  <  3q)  it  is  suggested  that  the  broadening  observed  on 
Fig.  C7  is  due  to  ARGUS  Event  2. 

Figure  C14,  a  3-dimensional  plot  of  true  count-rate  vs  B  and  L, 
summarizes  the  data  analyzed  for  the  unshielded  G-M  counter  (Detector  C, 
Channel  3)  for  ARGUS  Event  1. 

ARGUS  Event  2 

The  data  analyzed  for  ARGUS  Event  2  is  over  a  smaller  time  inter¬ 
val  than  Event  1.  Figure  C8  gives  the  true  count-rate  versus  elapsed 
time  since  Event  2.  Again  the  data  is  for  count-rates  above  the 
estimated  background.  The  peak  intensities  were  centered  veiy  nearly  at 
L  =  2.11.  Pass  453  was  the  first  to  indicate  ARGUS  Event  2,  0.03  hours 
after  the  event.  Considering  the  location  of  the  satellite  (297°  E, 

27°  N)  and  the  detonation  (352°  E,  50°  S),  fission  electrons  would  need 
to  have  energies  in  excess  of  19  MeV  to  reach  the  satellite  position  in 


their  eastward  drift  in  the  108  seconds  (Ref.  C25),  Although  there 
is  a  significant  numoer  of  MeV-range  8-electrons  available  from  the 
8-decay  spectra  of  fission  products,  the  8-decay  end  point  energies  are 
well  below  19  MeV.  Hence  the  initial  count-rate  may  be  due  to  direct 
illumination  by  'y  -rays  emitted  by  the  fission  debris.  The  bomb 
debris  would  have  to  expand  to  a  higher  altitude  than  the  detonation 
point  (about  480  km)  in  order  to  be  in  the  line  of  sight  of  the  satel¬ 
lite.  Colgate  (Ref.  C9)  discusses  the  mass  motion  of  a  high  altitude 
nuclear  explosion.  He  assumes  that  the  detonation  products  produce  a 
cavity  in  the  magnetic  field  and,  in  general,  the  hole  will  expand 
until  the  work  done  against  the  magnetic  field  equals  the  kinetic  energy 
of  the  detonation  products.  When  the  expansion  of  the  detonation  plasma 
products  ceases,  the  subsequent  motion  will  be  a  spreading  along  the 
lines  of  force  away  from  the  detonation  hemisphere.  This  occurs  in  the 
order  of  seconds  according  to  Colgate.  With  the  debris  distributed 
along  the  lines  of  force  for  L  =  2.11,  the  satellite  would  be  in  a 
position  to  detect  the  emitted  y  -rays. 

The  decay  rates  of  ARGUS  Event  2  are  very  nearly  represented  by 
a  1/t  law  if  the  data  for  Passes  466-492  are  considered.  However, 

Passes  505  and  518  indicate  increases  in  count-rate  in  a  manner  diffi¬ 
cult  to  explain  (Fig.  C9) .  A  sudden  commencement  of  solar  activity 
did  occur  shortly  after  Pass  505  and  3-5  September  are  listed  as  mag¬ 
netically  disturbed  days. 

It  is  noted  that  the  shielded  G-M  counter  (Detector  D,  Channel  1) 
had  very  low  count-rates  for  the  sequence  of  Passes  analyzed  for  ARGUS 
Event  2.  These  were  typically  less  than  20  see Figure  CIO  shows 
the  count-rate  ratio  of  the  unshielded/shielded  G-M  counter  for  the  L 
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shell  of  peak  intensity  (L  =  2.11)  for  Passes  466-492.  Subsequent  to 
Pass  492  the  ARGUS  Event  2  electrons  detected  by  the  shielded  counter 
were  indistinguishable  from  the  natural  radiation. 

C4.  ARGUS  Loss  Mechanisms 

The  effect  of  the  atmosphere  on  the  decay  of  the  trapped  electrons 
has  been  calculated  using  a  Fokker-Planck  sr>oroach  to  derive  a  diffusion 
equation  (Ref.  CIO).  This  technique  was  applied  by  :>7?lt  (Ref.  Cll),  also 
Walt  and  Newkirk  (Ref.  C12)  to  the  observed  time  variation  produced  by 
the  Starfish  nuclear  explosion.  In  regions  of  the  magnetosphere  for 
which  L  =  1.3,  it  was  found  that  the  dominant  factor  limiting  the 
trapping  time  of  electrons  is  scattering  by  components  of  the  atmosphere. 
However,  for  L  =  1.3  the  decay  of  trapped  electrons  is  dominated  by 
other  loss  mechanisms  of  electromagnetic-hydromagnetic  character. 

A  breakdown  of  the  adiabatic  invariants  which  will  eventually 
enhance  the  atmosphere  losses  has  been  proposed  by  a  number  of  authors. 
The  processes  responsible  for  the  violation  of  the  third  adiabatic 
invariant  have  been  considered  by  Parker  (Ref.  C13),  Davis  and  Chang 
(Ref.  C14) ,  and  Nakada  and  Mead  (Ref.  C15) .  Large  scale  geomagnetic 
fluctuations  associated  with  sudden  commencements  and  the  initial  phase 
of  geomagnetic  storms  can  cause  diffusion  of  electrons  across  L  shells. 
However,  the  dominant  effect  would  be  in  the  outer  radiation  zone, 
and  repeated  magnetic  storms  would  seem  to  be  necessary  for  rapid  dif¬ 
fusion.  The  exact  analytic  details  are  quite  complicated  and  have  not 
been  attempted  so  far  in  the  present  work  -  however  the  concept  is 
highly  suggestive.  Sudden  commencements  were  observed  on  August  27 
at  0243,  0303,  and  0541  hrs  UT.  Pass  414  traversed  the  ARGUS  Event  2 
region  at  approximately  0456  UT,  1.8  hours  after  the  event.  No  further 
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sudden  commencements  were  observed  until  September  3  and  4  which  co¬ 
incide  with  Passes  505  and  518  respectively. 

The  effect  of  several  other  electromagnetic- type  loss  contribu¬ 
ting  factors  have  been  studied  by  workers  in  this  field,  as  follows. 
Resonant  particle-wave  interaction  in  which  the  magnetic  moment  in¬ 
variance  is  violated  has  been  studied  by  Wentzel  (Ref.  C16),  Dragt 
(Ref.  C17),  and  also  Chang  and  Pearlstein  (Ref.  C18) .  A  summary  of 
Dragt ‘s  work  (Ref.  C-17)  on  the  effect  of  hydromagnetic  waves  is  given 
by  Fig.  Cll.  Critical  energies  below  which  the  magnetic  moment  is  con¬ 
served  in  interactions  with  hm  waves  are  shown  for  particles  mirroring 
at  latitudes  greater  than  30  degrees  assuming  the  hm  wave  frequency  is 
cut  off  sharply  at  v*  max»  Various  cut-off  frequencies  are  indicated 
as  it  is  assumed  that  hm  waves  may  vary  from  1  to  5  sec  *  up  to  the 
local  proton  cyclotron  frequency.  It  is  clear  that  such  frequencies 
would  not  be  very  effective  in  interaction  with  ARGUS  electrons. 

Kentzel  gives  a  qualitative  discussion  to  the  effect  that  hm  waves  could 
influence  the  mirror  point  distribution  of  the  relativistic  ARGUS 
electrons  during  a  geomagnetic  storm.  This  may  be  a  contributing 
factor  to  the  large  loss  observed  within  the  first  day  which  was  active 
magnetically. 

Again,  the  precipitation  of  electrons  through  the  action  of 
whistlers  has  been  investigated  by  Dungey  (Ref.  C19)  and  by  Cornwall 
(Ref.  C20) .  For  resonant  interaction  to  take  place,  the  Doppler- 
shifted  frequency  of  the  whistler  as  seen  by  the  electron  must  be  equal 
to  the  relativistic  cyclotron  frequency  of  the  moving  electron.  Figure 
C12  summarizes  Cornwall's  results  showing  the  L  dependence  of  the 
precipitating  particles  with  increasing  electron  energy.  The  ARGUS 
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Events  1  and  2  regions  a^e  shown,  indicating,  at  least  qualitatively, 
that  whistlers  ir.ay  not  be  a  major  mechanism  driving  electrons  to 
larger  B  values. 

Also,  Raulikas  et  al  (Ref.  C21)  suggest  longi'.udinal  waves 
propagating  along  magnetic  field  lines  would  interact  with  electrons  at 
their  mirror  points  and  would  be  an  efficient  loss  mechanism.  These 
electrostatic  waves  may  be  able  to  move  electrons  to  larger  B  values 
in  times  of  the  order  of  days. 

Further  detailed  theoretical  investigations  are  needed  before 
quantitative  results  may  be  given  to  adequately  explain  the  rapid 
decay  of  the  ARGUS  electron  shells. 

C5.  Conclusions 

The  analysis  of  the  data  obtained  by  the  omnidirectional  G-M 

counters  (Detectors  D,  C;  Channels  1,  3)  carried  by  Explorer  4  for  the 

ARGUS  Events  1  and  2  can  be  summarized  as  follows. 

The  two  Events  are  centered  on  the  L  =  1.72  and  L  =  2.11  shells 

respectively.  The  locations  of  the  peak  intensities  did  not  vary,  at 

least  within  experimental  error,  from  these  values  during  the  time  of 

-1  13 

observation.  The  ARGUS  Event  1  region  showed  a  decay  law  of  t 

for  E  >  3  MeV.  The  count-rate  ratio  of  the  unshielded  to  the  shielded 
e 

counter  was  found  to  decrease  with  time.  This  indicates  spectral 
changes  occur  which  is  in  apparent  contradiction  with  the  results  of 
the  published  Jason  rocket  data.  A  general  broadening  of  the  ARGUS 
Event  1  region  was  observed  subsequent  to  the  detonation  of  ARGUS 
Event  2.  The  increase  in  the  width  at  half-maximum  was  by  a  factor  of 
1.5.  ARGUS  Event  2  seemed  to  be  influenced  to  a  greater  degree 
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than  ARGUS  Event  1  by  the  magnetically  disturbed  days  of  September  3-S. 
Mcllwain  (Ref.  C22)  theorized  that  magnetic  storms  produce  a  decrease 
in  electron  flux  which  occurs  so  quickly  that  it  is  difficult  to  follow 
and  then  is  followed  by  an  increase  larger  than  normal.  The  observed 
increases  on  Passes  505  and  518  of  Event  2  may  be  related  to  these 
magnetic  storm  effects.  The  larger  value  L  shells  of  Event  2  compared 
to  Event  1  may  account  for  the  variation  of  the  storm  influence.  Al¬ 
though  a  general  decrease  in  count-rate  ratios  were  observed  for  both 
events,  ARGUS  Event  2  was  not  as  sharp  as  the  first  event  and  may  point 
to  L  dependence  of  the  decaying  fission  spectrum.  The  data  is  too 
limited  to  determine  the  precise  nature  of  this  dependence. 

The  decay  of  the  ARGUS  Event  1  region  was  found  to  follow  a 

-1  13 

t  '  law.  It  may  be  coincidental  that  this  is  nearly  identical  to 

the  fission  product  decay  law.  The  fission  fragment  p-ray  source 

-1  2 

strength  decays  approximately  according  to  a  t  '  law  (Ref.  C23). 

-1  2 

Measurements  indicate  that  although  the  t  *  decay  represents  a 

reasonable  average  there  have  been  instances  where  exponents  in  the 

range  -0.9  to  -2.0  are  required  to  represent  the  rate  of  decay. 

Figure  13,  from  Heller  (Ref.  C24),  gives  the  number  of  beta  rays  per 
235 

fission  from  U  per  energy  interval  versus  time  for  various  energies. 
The  ratio  of  high  energy  to  low  energy  y3-rays  decreases  with  time. 

This  contradicts  the  results  of  Figs.  C6  and  CIO  which  indicated  the 
unshielded/shielded  count-rate  ratio  that  is,  C(Ee  >3  MeV)/ 

C(Eg  >  7.6  MeV).  decreases  with  time.  Further  study  is  required  on 
the  lifetimes  of  the  fission  debris  which  is  trapped  in  the  earth's 
geomagnetic  field. 

A  brief  study  was  conducted  of  possible  mechanisms  which  would 
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explain  the  rapid  decay  of  the  ARGUS  Events.  Theories  involving  the 
breakdown  of  the  adiabatic  invariants  have  been  proposed  by  a  number 
of  workers.  The  author  (J.  George)  applied  these  theories  to  the 
ARGUS  data.  Computations  indicate  that  hydromagnetic  waves  and  whistler 
interactions  resulting  in  the  violation  of  the  first  invariant  are  not, 
at  least  qualitatively,  major  mechanisms  in  removing  ARGUS  electrons. 
Although  resonant  particle-wave  interactions,  and  other  processes  such 
as  diffusion,  may  operate,  they  do  not  dominate.  Further  theoretical 
work  is  necessary  to  determine  the  controlling  mechanism  for  the  ARGUS 
decay. 

Adiabatic  processes  such  as  betatron  acceleration  and  decelera¬ 
tion  cannot  be  excluded  in  any  detailed  analysis  of  irreversible  loss 
mechanisms.  The  effect  of  slow  changes  in  the  magnetic  field  upon  the 
omnidirectional  intensity  of  particles  was  considered  and  calculations 
for  the  ARGUS  time  period  show  that  the  changes  are  not  inconsequential. 
Ry  making  conservative  assumptions  concerning  the  change  in  the  magnetic 
field  during  the  main  phase  decrease  following  Pass  414,  it  was  shown 
that  the  relative  intensity  could  change  by  15%  due  to  the  betatron 
acceleration  alone. 
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Table  Cl  (Ref.  C2)  Characteristics  of  G-M  counters,  Explorer  4 


Detector 


Shielding 


Sensitive 

to: 


Geometric 

Factor: 


D  (Channel  1) 

"Shielded1' 
Geiger-Mueller 
counter  (Anton 
302)  Cylinder 
approximately 
7  mm  x  9  mm 
Scalar:  64 


1.2  gVcm  Fe 
+1.6  gm/cm2 
Pb  (minimum) 


Electrons  of 
Ee>  5  MeVj 
Protons  of 
E^40  MeV, 
X-rays  of 
E>  80  KeV 
with  low 
efficiency 


Omnidirectional 
geometric  factor  G  : 
0.14  cm2  for  0 
minimum  stopping 
power,  0.82  cm2 
for  penetrability 
greater  than  7 
gm/cm2 


C  (Channel  3) 

"Unshielded" 
Geiger-Mueller 
counter  (Anton 
302)  Cylinder 
approximately 
7  mm  x  9  nun 
Scalar:  2048 


2 

1.2  gm/cm  Fe 
(minimum) 


Electrons  of 
\>  3  MeVj 
Protons  of 
E l>  30  MeV, 
X-rays  of 

E>20  KeV 
with  low 

efficiency 


Omnidirectional 
geometric  factor  G  : 
0.14  cm2  for  0 
minimum  stopping 
power,  0.70  cm2 
for  penetrability 
greater  than  5 
gm/cm2 
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Table  C2.  Shielding  and  Range-energy  data  for  G-M  Detector  C,  D 

Detector 

Shielding 

Equivalent  Al. 

Energy 

(gVcm2) 

2 

(gVcm  ) 

(MeV) 

C 

1.2  Fe 

1.46 

2.95 

D 

1.2  Fe 

+  1.6  Pb 

5.89 

| 

7.6 
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Table  C3.  Explorer  4  Passes  ''over'-  Huntsville,  Alabama  selected  for 
omnidirectional  radiation  analysis 


ARGUS  Event  1  (time  of  burst:  0230  UT  August  27,  1958) 


Date 

Pass 

Pass 

Pass 

August  27 ,  1958 

414 

415 

4 16 

28 

427 

29  • 

440 

30 

453 

454 

31 

466 

September  1,  1958 

479 

2 

492 

3 

505 

4 

518 

ARGUS  Event  2  (time  of  burst:  0320  UT  August  30,  1958) 


Date 

Pass 

Pass 

Pass 

August  30,  1958 

453 

454 

455 

31 

466 

September  1,  1958 

479 

2 

492 

3 

505 

4 

518 

Cl6 


2 

Pig.  Cl.  Geometrical  factor,  Gq  (cm  ),  versus  the  energy  (MeV) 
necessary  for  transmission  by  the  various  thicknesses 
of  materials.  Circled  points  are  from  Table  C2,  using 
Eq.  (Cl),  assuming  a  linear  dependence  of  Gq  with  E 
between  plotted  points  and  constant  values  for  larger 
energies  (as  shown). 
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Pig.  C2.  Omnidirectional  count-rates  vs  time  for  a  typical  Pass  of 

Explorer  4  in  the  ARGUS  Event  1  region.  The  solid  lines 

are  average  count-rates  drawn  through  points  not  shown. 

Curve  A.  Unshielded  G-M  counter  (Detector  C,  Channel  3) 
which  detects  electrons  of  energy  greater  than 
3  MeV.  The  background  has  been  drawn  in  by 
estimation. 

Curve  B.  Shielded  G-M  counter  (Detector  D,  Channel  1  ) 

which  detects  electrons  of  energy  greater  than 
7.6  MeV.  The  background  has  been  drawn  in  by 
estimation. 

Curve  C.  Ratio  of  A/B,  unshielded  to  shielded  G-M  count- 
rates.  The  ratios  are  of  count-rates  above  the 
background . 

Curve  D.  Difference  A  -  B,  unshielded  minus  shielded  G-M 
count-rates  above  the  background.  This  curve 
shows  the  omnidirectional  count -rates  due  to 
electrons  between  3  and  7*6  MeV,  a  roughly 
"monoenergetic"  count-rate  curve. 
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Pig.  C3. 


Patn  of  satellite  Explorer  4  in  B-L  coordinates, 
regions  of  ARGUS  Events  1  and  2  are  shown  also. 


The 
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Fig.  C4, 


Planetary  magnetic  three-hour 


range  K 
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indices  during 


the  time  interval  of  ARGUS  Events  1  and  2. 
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Fig.  C$.  Logarithmic  plot  of  true  count-rate  of  the  unshielded 
G-M  counter  (Detector  C,  Channel  3)  versus  the  elapsed 
time  since  detonation  of  ARGUS  Event  1.  The  solid  lines 
are  least  squares  fit  of  Eq.  (C3)  for  the  data  of  the 
sequence  of  Passes  414,  427,  •  ••  5*8 .  n  =  1.13  for 
both  curves.  Intersection  of  the  two  straight  lines 
with  the  ordinate  axis  yield  values  of  CQ,  Eq.  (C3). 

The  peak  intensities  were  centered  at  L  =  1.72. 
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Pig.  C6.  Count-rate  ratio,  unshielded  {Channel  3)/shielded 

(Channel  1)  omnidirectj.onal  G-M  counters  versus 

elapsed  time  Te  since  ARGUS  Event  1.  Both  G-H 

count-rates  were  corrected  for  background  count- 

rates.  The  dashed  line  is  the  plot  of  C  =  C  t 

0 

Eq.  (C3). 
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RATIO 


1=1.72 


Pig.  C7.  The  geometrical  width  at  half -maximum  of  the  ARGUS 
Event  1  region  versus  the  elapsed  time  Tg  since  the 
event.  (Pass  453  was  the  first  to  contain  ARGUS 
Event  2  data.) 
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Fig.  C8.  True  count-rate  of  the  unshielded  (Detector  G, 

Channel  3)  G-M  counter  versus  the  elapsed  time 

T  since  ARGUS  Event  2  for  selected  L  shells, 
e 
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Pig.  C9.  Logarithmic  plot  of  true  count-rate  of  the  unshielded 

G-M  counter  (Detector  C,  Channel  3)  vs  the  elapsed 

time  T  since  ARGUS  Event  2  for  selected  L  shells, 
e 

Passes  505  and  518  were  during  a  period  of  high  magnetic 
activity. 
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Pig.  CIO.  Count-rate  ratio,  unshielded  (Channel  3)/shielded 
(Channel  1)  G-M  counter  vs  elapsed  time  Tq  since 
ARGUS  Event  2.  L  =  2.11  is  the  region  of  peak 
intensity  of  the  ARGUS  Event  2  region.  The  numbers 
above  the  plotted  points  are  satellite  Pass  (over 
Huntsville,  Alabama)  numbers. 
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Pig.  Cll.  Critical  energies  below  which  the  magnetic  moment 

is  conserved  in  interactions  with  hydromagnetic 

waves  for  particles  mirroring  at  geomagnetic 

latitudes  greater  than  30°  assuming  the  hm  wave 

is  cut  off  sharply  at  ^max.  Shown  are  Hnax  = 

1,  5,  and  V,  sec"1  (Ref.  C19). 

1* 


C2? 


^critical  (MeV> 


L  (earth  radii) 


Pig.  C12.  Regions  in  space  where  whistlers  may  be  effective 
in  precipitating  electrons  into  the  atmosphere. 
The  plot  is  from  Eq.  (C4)  with  the  heart  of  the. 
whistler  range  being  1  to  10  kc/s  (Ref.  C20). 
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Pig.  C13.  Number  of  beta  rays  per  fission  from  U  per  energy 

? 

interval  for  various  energies  E  in  units  of  mQc  . 

The  data  is  taken  from  Heller  (Ref.  C2<f). 
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Fig.  Cl4.  Observed  omnidirectional,  true  initial  count-rate  CQ 

contours  vs  B  and  L  for  ARGUS  Event  1  (Detector  C, 

unshielded  G-M  counter,  Channel  3).  C  is  the  true 

o 

-1  13 

count-rate  corrected  for  decay  assuming  the  t 
law,  Eq.  (C3).  The  Initial  time  reference  for  the 
decay  law  correction  is  the  beginning  of  ARGUS  shell 
counting  for  Pass  427.  The  solid  line  count-rate 
contour  curves  are  best  fits  to  the  data  from  Passes 
415,  4l6,  427,  440,  453,  466,  479  (points  not  shown). 
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Pig.  C15.  Omnidirectional  flux  density  JQ  contours  vs  B 

and  L  for  ARGUS  Event  1  (Detector  C,  unshielded 
G-H  counter,  Channel  3)  obtained  from  the  observed 
omnidirectional  true  count-rate  CQ  data  (Pig.  Cl4). 
The  conversion  is  given  by  Jq  =  67. Cq. 
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D.  Summary  of  M.S.  Thesis,  Joseph  R.  Ockersz,  S.J.,  "Omnidirectional 
Flux  Densities  of  Geomagnetically  Trapped  Radiation  Explorer  4 
Pre-ARGUS  Times";  Saint  Louis  University,  Physics  Department 
(May  1966). 

Dl.  Introduction 

The  purpose  of  the  present  Summary  is  to  present  experimental 
aspects  of  the  trapped  radiation,  including  the  omnidirectional  charged 
particle  flux  densities  in  the  geomagnetic  field,  for  pre-ARGUS  event 
times  of  Explorer  4  Satellite  data  telemetered  over  Huntsville,  Alabama. 
Such  an  investigation  is  important  toward  understanding  the  "background" 
radiation  (natural  radiation)  against  which  the  artificially  injected 
ARGUS  radiation  is  observed. 

D2.  Experimental  Details 

Instrumentation  and  Calibration  of  Explorer  4 
The  University  of  Iowa  group  designed  the  detection  instruments. 
Van  Allen  has  published  a  complete  description  of  the  detection  instru¬ 
ments  (Ref.  D2) ,  Reference  D3,  Fig.  2,1,  p.  16,  shows  a  block  diagram 
of  the  four  detectors  A,  B,  C  and  D,  and  the  corresponding  telemetry 
channels  2  and  5,  3  and  1.  (See  also  Preface  P2,  Fig.  P2-1). 

For  the  present  thesis  the  omnidirectional  count-rates  were 
analysed,  involving  data  obtained  by  G-M  Detectors  C  and  D.  Preface 
P2  of  this  Report  summarizes  the  characteristics  of  Detectors  C  and  D. 
Table  Dl  lists  the  telemetry  data  available  for  the  present  analysis 
and  Fig.  Dl  shows  the  location  of  Satellite  Explorer  4  for  some  of  the 
Passes . 

The  count-rates  of  the  G-M  counters  must  be  corrected  for 


Dl 


"deadtimes".  For  detectors  C  and  D,  Anton  502  G-M  co-inters,  the  cqua 
tion  for  the  correction  is 


C  =  Ce^C 
obs 


(Dl) 


where  C  ,  is  the  observed  count-rate  and  C  is  the  true  or  corrected 
obs 

count-rate  and  7«  ,  the  deadtime,  is  62.5_+1.3  microseconds. 

Further  instrumentation  and  analysis  details  are  found  in  DASA 
Report  WT  1613,  Ref.  (D3) . 


D3.  Results 

A  complete  set  of  Figures  showing  the  true  count-rates  vs  time 
for  the  different  B-L  regions  of  the  pre-ARGUS  Explorer  4  times  is 
found  in  Ref.  Di. 

These  Ref.  Dl  Figures  are  here  compactly  summarized  by  three  new 
Figs.  D2  (constant  omnidirectional  count-rate  contours  in  B-L  coordi¬ 
nates)  ,  D3  (omnidirectional  flux  densities  vs  B) ,  D4  (three  dimensional 
plots  of  omnidirectional  true  count-rates  vs  B  and  L) . 

The  slope  of  Log(«jT)  vs  Log(B)  straight  line  (derived  from 
Fig.  D4)  compares  closely  with  earlier  published  results  (Ref.  D4). 

In  general  it  has  been  found  (Ref.  Dl)  that  Passes  2,  13,  26,  52,  78, 
107,  173,  207  and  367  recorded  small  flux  densities,  between  10  and  30 
particles  cm'^sec"*;  while  passes  4,  43,  121,  272,  285,  298,  311,  324, 

337,  375  and  401  recorded  large  flux  densities  between  100  and  20,000 
-2  -1 

particles  cm  sec  .  Passes  of  small  measured  flux  densities  correspond 
to  low  altitude  (250  <  altitude  <  700  km) .  Tne  low  count-rates  are 
due  chiefly  to  non- trapped  cosmic  radiation  in  space.  The  Passes  of 


D2 


larger  measured  altitudes  (1000  altitude  <  1500)  are  specified  by 

0. 17  <  B  <  C.30,  i.3S  <  L  <  1.65  and  0.35  <  B  <0.45,  3.5  <  L  ^  5.5. 

The  larger  omnidirectional  flux  densities  were  recorded  in  the  Van 

Allen  region. 
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Fig.  D1  Map  shov/ing  location  of  Satellite  1958€(Explorer  4)  for 

Passes  2,  4,  13,  26,  43,  52,  63,  78,  107,  121,  172,  207,  272, 
285,  298,  311,  324,  337,  367,  375,  401  and  414  tracked  by  the 
U.  S.  Army  Tracking  Station  at  Huntsville,  Alabama.  Except 
for  Pass  414,  these  Passes  are  for  pre-ARGUS  times.  The 
Turber  in  parenthesis  is  the  Pass  number,  the  other  numbers 
are  the  satellite  altitudes  at  the  beginning  and  end  of  the 
tracking  period  for  the  Pass. 


Fig.  D2  Omnidirectional  true  constant  count-rate  contours  in  B-L 
coordinates.  Data  is  for  pre-A'IGUS  times,  shielded  G-M 
detector  D  (channel  1)  for  the  region  defined  by 
1.27  <  L  <  1.30  and  0.20  <  B  <  0.23.  To  obtain  Jq,  the  true 
omnidirectional  count-rates  are  multiplied  by  8.4*,  a  reduc¬ 
tion  constant  derived  from  the  calibration  of  the  detector, 
(♦provisional,  subject  to  correction) 
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Fig.  D3  Omnidirectional  flux  densities  vs  magnetic  flux  density  B. 

Data  is  for  pre-ARGUS  times,  shielded  G-M  detector  D  (channel  1) 
for  the  region  defined  by  0.20  <■  B  <  0.23  and  1.27  <  L  <  1.30. 
Log  JQ  vs  B  is  plotted  instead  of  log  JQ  vs  log  B  in  order  to 
obtain  somewhat  greater  spread  in  data.  The  data  may  be  repre¬ 
sented  by  the  equation  =  k  Bn,  whence  n  =  30,  approximately. 
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Fig.  D4  Three  dimensional  plot  of  omnidirectional  true  count-rates  C 
(shielded  G-M  counter.  Detector  D,  channel  1)  (number  sec"*) 
vs  B  and  L.  Sheets  of  constant  L  are  shown  intersecting  the 
contours  of  constant  count-rates.  Data  is  the  same  as  for 
Figs.  D2  and  D3. 


Table  Dl.  Explorer  4  telemetry,  pre-ARGUS  times  data  available  at 
Saint  Louis  University,  Physics  Department;  Fasses  are 
over  Huntsville,  Alabama  Tracking  Station. 


Pass 

No. 

Date 

Time 

1 

July  -  26 

1658-1702 

2 

26 

1855-1858 

3 

26 

2054-2057 

4 

26 

2240-2257 

5 

July  -  27 

0036-0042  ! 

13 

27 

1650-1655 

15 

27 

1348-1851 

16 

27 

2042-2052 

17 

27 

2240-2247 

26 

July  -  28 

1449-1451 

27 

28 

1642-1648 

28 

28 

1840-1844 

42 

July  -  29 

2026-2035 

43 

29 

2220-2227 

52 

July  -  30 

1430-1437 

53 

30 

1627-1632 

64 

July  -  31 

1421-1426 

66 

31 

i6l8-l622 

67 

31 

1816-1818 

69 

31 

2205- 

70 

Aug.  -  1 

0000-0007 

77 

1 

1220-122 2 

78 

1 

1412-1420 

79 

1 

1609-1614 

80 

1 

1806-1811 

81 

1 

2001-2008 

82 

1 

2154-2208 

94 

Aug.  -  2 

1951-1957 

104 

Aug.  -  3 

-0841 

107 

3 

1941-1949 

108 

3 

2130-2141 

109 

3 

2331-2335 

121 

Aug.  -  4 

2125-2134 

122 

4 

2321-232? 

133 

Aug.  -  5 

1921-1927 

134 

5 

2113-2126 

147 

Aug.  -  6 

2200-2210 

148 

6 

1257-1305 

I69 

Aug,  -  ? 

1257-1305 

172 

Aug.  -  8 

1845-1852 

173 

8 

2039-2050 

182 

Aug.  -  9 

1243-1249 

189 

9 

1820-1828 

Pass 

No. 

Date 

Time 

199 

Aug.  -  10 

2015-2025 

207 

11 

1025-1035 

208 

11 

1220-1228 

220 

Aug.  -  12 

1010-1021 

221 

12 

1807-1815 

222 

12 

1405- 

22  + 

12 

1754-1801 

225 

12 

1948-1957 

233 

Aug.  -  13 

0957-1008 

234 

13 

1153-1202 

237 

13 

1741-1748 

238 

13 

1935-1943 

246 

Aug.  -  14 

0950-0953 

24? 

14 

1139-1149  ! 

250 

14 

1725-1734 

251 

14 

1923-1929 

259 

Aug.  -  15 

0924-0938 

260 

15 

1125-1132 

264 

15 

1907-1915 

2?  2 

Aug.  -  16 

0910-0920 

273 

16 

1107-1118 

276 

16 

1657-1705 

277 

16 

I852-I858 

285 

Aug.  -  17 

0855-0906 

286 

17 

1056-1102 

298 

Aug.  -  18 

0840-0854 

302 

18 

1627-1633 

310 

Aug.  -  19 

0630-0642 

311 

19 

0823-0833 

312 

19 

1022-1032 

313 

19 

I6ll-l6l7 

324 

Aug.  -  20 

0836-0809 

325 

20 

1004-1015 

328 

20 

1555-1602 

337 

Aug.  -  21 

0748-0806 

338 

21 

0949-0959 

341 

21 

1538-1544 

355 

Aug.  -  22 

1521-1527 

367 

Aug.  -  23 

1503-1510 

375 

Aug.  -  24 

0510-0518 

401 

Aug.  -  26 

0433-0440 

406 

26 

1408-1414 

413 

26 

2225-2259 
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E.  Summary  of  Pertinent  Analytic  Sections  of  Ph.D.  Dissertation, 

Joseph  M.  Palkeday,  "Interpretation  of  Directional  Flux  Densities 
in  ARGUS  Shells,  Explorer  4  Satellite  Data";  Saint  Louis  Univer¬ 
sity,  Physics  Department  (March  1966). 

El.  Introduction 

The  method  used  for  obtaining  the  directional  flux  density  J ( 0 1 ) 
-2  -1  -1 

(particles  cm  sec  ster  ;  G!,  angle  between  scintillation  counter 

axis  P  and  J)  follows  that  of  Lundquist,  Nauraann,  and  Weber  (Ref.  El) 

in  the  reduction  of  the  integral  equation  relating  the  corrected  direc- 

-2  -1 

tional  count-rate  C ( 0)  (number  cm  sec  )  and  J(0')  to  matrix  form 

A 

(0,  angle  between  P  and  a  plane  perpendicular  to  the  magnetic  flux 
density  vector  B).  The  matrix  elements  S^  (Ref.  E2,  p.  13)  associated 
with  the  response  function  have  been  evaluated  to  a  greater  accuracy 
resulting  from  a  more  precise  definition  in  the  present  analysis  of  the 
limits  of  Integration  in  the  Integral  equation  relating  C(0)  and  J(0‘) 
(Ref.  E2,  Eq.  (1.6)).  It  has  been  found  that  the  kernel  of  this  inte¬ 
gral  equation  has  a  vanishingly  small,  determinant  and  as  a  consequence 
established  methods  of  solving  the  integral  equation  for  J ( 0* )  values 
fail  to  produce  physically  meaningful  results.  The  kernel  of  this 
integral  equation  is  highly  sensitive  to  experimental  errors  inherent 
in  the  C ( 0 )  data.  The  present  method  replaces  the  inverse  of  the 
kernel  matrix  by  an  approximate  matrix  (Ref.  E2,  Eq.  (1.26))  which  is 
insensitive  to  small  errors  in  the  counting  rate  data. 


El 


E2.  Summary  of  the  Analytic  Procedure 

The  complete  details  of  the  analytic  procedure  outlined  in  El 
above  are  given  in  the  Paikeday  Ph.D.  Dissertation  (Ref.  E2/.  The 
procedure  involves  the  following  steps: 

(1)  Reduction  of  the  integral  equation  to  a  matrix  form  (see  Ref.  E2 
PP.  7-13). 

(2)  Representation  of  the  directional  scintillation  detector  (Detector 
A,  Channel  2)  response  function  S(a)  (Ref.  E2,  Pig.  1.2,  p.  149) 
in  analytic  form  (see  Ref.  E2,  Eq.  (1.10))  to  yield  J ( 0 * )  values. 

(3)  Inversion  of  the  matrix  equation  (see  Ref.  E2,  pp.  14-33). 

E3.  Conclusions 

With  the  data  available  as  discrete  points  of  the  C(Q)  function, 
the  solutions  obtained  by  using  the  inverse  of  the  S-matrix  are  found 
to  be  non-physical.  The  process  of  getting  correct  physical  solutions 
does  not  in  any  way  remove  the  ill-conditioning  of  the  original  matrix 
S  representing  the  kernel  of  the  integral  equation  but  rather  supplies 
information  about  the  type  of  errors  that  enter  the  counting  rate  data. 

The  method  developed  in  the  present  analysis  is  found  to  yield 
meaningful  J ( 0 ' )  values  for  the  available  C ( 0)  data. 
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